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A massive black hole is present at the centre of our galaxy and inevitably accretes 
dark matter particles, creating a region of very high particle density. The annihi- 
lation rate is enhanced with a large number of e+e~ pairs produced either directly 
or by sucessive decays of mesons. We evaluate the synchrotron emission (and self- 
absorption) associated with the propagation of these particles through the galactic 
magnetic field and give constraints on the values of mass and cross section of the 
dark matter particles. 

1 Introduction 

There is convincing evidence for the existence of an unseen non-baryonic com- 
ponent in the energy density of the universe. The most promising dark matter 
candidate appear to be weakly interacting massive particles (WIMPs) and in 
particular the so-called neutralinos, arising in supersymmetric scenarios (for 
a review of particle candidates for dark matter see e.g. El). The annihilation 
of these X-particlcs would produce quarks, leptons, gauge and Higgs bosons 
and gluons. In particular e^e" pairs are produced either directly or by suc- 
cessive decays of mesons, and lose their energy through synchrotron radiation 
as they propagate in the galactic magnetic field. This radiation is expected 
to be greatly enhanced in the proximity of the galactic centre, where the exis- 
tence of a massive black hole creates a region of very high dark matter particle 
density and consequently a substantial increase in the annihilation rate and 
in the ensuing synchrotron radiation. 

2 Dark matter distribution around the Galactic Centre 

There is strong evidence for the existence of a massive compact object lying 
within the inner 0.015pc of the galactic centre (see □ and references therein). 
This object is a compelling candidate for a massive black hole, with mass M = 
2.6±O.2M0. The galactic halo density profile is modified in the neighbourhood 
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Figure 1. Left panel: A,^ as a function of frequency for mx = ITeV. The two upper curves 
correspond to the cross section crv = 10~^*/m,^ (GeV) cm'^s"^, close to the unitarity limit; 
the two lower curves correspond to av = 10~^*/m,^ (GeV) cm"^s~^, a cross section more 
typical for wimps. Results for two values of the density profile are shown in each case: 
7 = 1, for solid curves and 7 = 1.5 for dashed ones. Right panel: Au as a function of the 
particle mass for i^=408MHz, av = 10~'^*/r7i^ (GeV) cm'^s"-'^ and two values of 7. 



of the galactic centre by the adiabatic process of accretion towards the central 
black hole. If we consider an initially power-law type profile of index 7, as 
predicted by high resolution N-body simulations u, the corresponding dark 
matter profile after accretion is modified to 
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where = (9 — 27)/(4 — 7), D is the solar distance from the galactic centre 
and PD = Q.2A:GeV / / cm? is the corresponding density. The factors a-y and 
g-y(r) cannot be detemiined analytically (for approximate expressions and 
numerical values see □). Expression |l] is only valid in a central region of 
size Rsp — ajD{M/ pDD^y^?~''^ where the central black hole dominates the 
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Figure 2. Comparison of Sgr A* (see Narayan et ai(1998)) observed spectrum with expected 
fluxes. The values of particle mass and cross section were choosen to flt the experimental 
data normalisation. 



gravitational potential. If we take into account the annihilation of dark matter 
particles, the density cannot grow to arbitrarily high values, the maximum 
density being set by the value 

Pcore = m/avtBH (2) 

where Ibh ~ 10^°years is the age of the central black hole. The final pro- 
file, resulting from the adiabatic accretion of annihilating dark matter onto a 
massive black hole, is 

Pdm{r) = -TT— (3) 

P(r) + Pcore 

following a power-law profile for large values of r, and with a flat core of 
density pcore and dimension 



-^core — ^sp I ) (4) 
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3 Constraints from Synchrotron Emission 



Among the products of annihilation of dark matter particles, there are en- 
ergetic electrons and positrons, which are expected to produce synchrotron 
radiation in the magnetic field around the galactic centre. The pair produc- 
tion spectrum is determined by the quark fragmentation function and has 
been evaluated using the MLLA approximation (see eI for details) . The galac- 
tic magnetic field can be modeled by making the 'equipartition assumption', 
where the magnetic, kinetic and gravitational energy of the matter accreting 
on the central black hole are in approximate equipartition (see eI). In this 
case, the magnetic field can be expressed as 

Bir)^l,G[^y^' (5) 

Most of the annihilations occur at very small distances from the centre, 
typically at « min{Rcore, 10i?s), i.e. in a region with magnetic fields of the 
order of > IG. Under these conditions, comparable to the size of the region 
where most of the annihilations occur, the electrons lose their energy almost 
in place. To compute the synchrotron luminosity resulting from the propa- 
gation of in the galactic magnetic field, we need to evaluate their energy 
distribution in the magnetic field, following u 

dn rY,{> E) 

dE = -l^^^^'^ 

where T is the annihilation rate 

r = ^ / p%4nr' dr, (7) 
mi Jo 

the function fe{r) is given by 

/'^(^) = foo /7 , , (8) 
Jo Pip47rr^ dr 



and 

(9) 

3m|c' 

The quantity Ye{> E) is the number of e"*"e~ with energy above E produced 
per annihilation, which depends on the annihilation modes. The energy de- 
pendence of Y(.{> E) can be neglected, and we estimate Ye(> E) by the 
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number of charged particles produced in quark fragmentation (see u). We 
obtain a simple analytical expression for the total synchrotron luminosity 

1/2 
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where 



The synchrotron self- absorption coefhcent is defined 13 by 
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where t(6) is the optical depth as a function of the cylindrical coordinate b 
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and the coefhcent is given by 
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The final luminosity is obtained by multiplying eq. ( pO| ) with given by 
eq. (p^). It is evident that in the limit of small optical depths the coefficent 
Ai, — > 1, as can be seen by expanding the exponential. We find for a^, the 
following expression 

FY c2 

"'^ = 1 3 

47r 

which can in turn be used to evaluate r(&) in eq. ( [T^ ) and in eq. (^2|). 

We first evaluate the self-absorption coefficent for selected values of the 
mass, as a function of frequency. The coefficent (see left part of fig. |l|) grows 
from very low values and then reaches the value 1, around a frequency which 
is strongly dependent on the cross section and the mass mx, but not very 
much on the profile power-law index 7. The right part of fig.0 shows the self- 
absorption coefficent at the fixed frequency of408MHz as a function of the 
wimpzilla mass. The behaviour shown is qualitatively the same for any value 
of the cross-section and for different 7. 

The coefficent is evaluated for two different values of the cross sec- 
tion, the first one corresponding to the maximum possible value av » 
10~^®/m^(GeV) cni-^s"^, (the so-called unitarity bound, see E3) and the sec- 
ond one corresponding to typical cross sections for supersymmetric scenarios. 



proof: submitted to World Scientific on February 1. 2008 5 




>8 



Mx (GeV) 



Figure 3. Exclusion plot based on the comparison between predicted flux and radio obser- 
vations of the galactic centre. The 3 solid curves indicate, for 3 different density profile power 
law index, the lower edge of the excluded regions. The dashed line shows, for comparison, 
the unitarity bound, av ~ The shaded region is the portion of the parameter space 

occupied by cosmologically interesting neutralinos (i.e. those leading to 0.025 < ^Ixh^ < 1; 
see, e.g. Bergstrom et al (f997)). 



(jv = 10~'^*/TO^(GeV) cm'^s"^ (see, e.g., In figure ^ we compare the pre- 
dicted spectrum with the observations; we choose a set of parameters mx, 
7 and av in order to reproduce the observed normahsation. It is remarkable 
that in this way one can reproduce the observed spectrum over a significant 
range of frequencies. The submillimeter excess is generally attributed to pro- 
cessesJji the accretion disk within several Schwarzschild radii of the black hole 
(e.g. EJ). The set of dark matter parameters for which the fluxes predicted in 
our model are consistent with observations is shown in the exclusion plot of 
fig. I 
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4 Conclusions and perspectives 

We have shown that present data on the emission from Sgr A* are compatible 
with a wide set of dark matter parameters. The evaluation of synchrotron 
self-absorption has enabled us to reach an alternative conclusion to an earlier 
claim of incompatibility of cuspy halos with the existence of annihilating wimp 
dark matter. We find that the experimental data on the Sgr A* spectrum at 
radio wavelengths could be explained by synchrotron emission of electrons 
produced in the annihilation of rather massive dark matter particles, extend- 
ing up to and beyond TeV masses. This is relevant to a recent study of 
coannihilations t3, which suggests that WIMPs with flxh'^ = 0.2 can extend 
up to several TeV, as well as to very massive particles (wimpzillas) that are 
produced non-thermally in the primordial universe. We also find that the 
synchrotron emission tends to give somewhat sharper constraints on masses 
and cross sections than the observed gamma-ray fluxes and neutrino limits. 
This situation could change with more sensitive gamma-ray and neutrino ob- 
servations anticipated from forthcoming experiments. 
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